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Abstract 

We study hypoelliptic operators with polynomially bounded coefficients that are of 
the form K — + Xq + /, where the Xj denote first order differential 

operators, / is a function with at most polynomial growth, and Xf denotes the formal 
adjoint of Xi in \3 . For any e > we show that an inequality of the form < 
C(||M||o,e + \\{K + iy)it||o,o) holds for suitable 5 and C which are independent of 
J/ £ R, in weighted Sobolev spaces (the first index is the derivative, and the second the 
growth). We apply this result to the Fokker-Planck operator for an anharmonic chain 
of oscillators coupled to two heat baths. Using a method of Herau and Nier [HN02], 
we conclude that its spectrum lies in a cusp {x + iy\x > — c, r G (0, 1], c G R}. 

1 Introduction 

In an interesting paper, [HN02], Herau and Nier studied the Fokker-Planck equation 
associated to a Hamiltonian system H in contact with a heat reservoir at inverse tem- 
perature (3. For this problem, it is well-known that the Gibbs measure 



is the only invariant measure for the system. In their study of convergence under the 
flow of any measure to the invariant measure, they were led to study spectral proper- 
ties of the Fokker-Planck operator C when considered as an operator on l^ifip). In 
particular, they showed that C has a compact resolvent and that its spectrum is located 
in a cusp-shaped region, as depicted in Figure 1 below, improving (for a special case) 
earlier results obtained by Rey-Bellet and Thomas [RBT02b], who showed that e^'''* is 
compact and that C has spectrum only in ReA > c > aside from a simple eigenvalue 
atO. 

Extending the methods of [HN02], we show in this paper that the cusp-shape of the 
spectrum of C occurs for many Hormander-type operators of the form 



Hpidpdq) = exp {~(3H(p, q)) dpdq 




(1.1) 
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Figure 1 Cusp containing the spectrum of £. 

(the symbol ^ denotes the formal adjoint in L'^) when the family of vector fields 
{-^jVjLa is sufficiently non-degenerate (see Definition 2.1 and assumption bi below) 
and some growth condition on / holds. 

The main motivation for our paper comes from the study of the model of heat 
conduction proposed in [EPR99a] and further studied in [EPR99b, EHOO, RBTOO, 
RBT02b, RBT02a]. These papers deal with Hamiltonian anharmonic chains of point- 
like particles with nearest-neighbor interactions whose ends are coupled to heat reser- 
voirs modeled by linear classical field theories. Our results improve the detailed knowl- 
edge about the spectrum of the generator C of the associated Markov process, see 
Sect. 5. As a by-product, our paper also gives a more elegant analytic proof of the 
results obtained in [EHOO]. A short probabilistic proof has already been obtained in 
[RBT02b]. 

The main technical result needed to establish the cusp-form of the spectrum is the 
Sobolev estimate Theorem 4. 1 which seems to be new. 
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2 Setup and Notations 

We will derive lower bounds for hypoelliptic operators with polynomially bounded 
coefficients that are of the form (1.1). We start by defining the class of functions and 
vector fields we consider. 

2.1 Notations 

For iV e R, we define the set Pol^ of polynomially growing functions by 

PoC = 1/ e C°°(R") Va, sup (1 + \\x\\)-^d"fix) < cA . (2.1) 
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In this expression, a denotes a multi-index of arbitrary order We also define the set 
Poll of vector fields in R" that can be written as 

n 

G = Goix) + J2 Gj(x)dj , Gi e Pol^ . 
j=i 

One can similarly define sets Pol^ of fcth order differential operators. It is clear that if 
X G Pol^ and Y G Polf , then [X, Y\ G Polf_+^i^i ■ If / is in Pol^, but not in Pol^+^ 

for any e > 0, we say it is of degree iV. 

2.2 Hypotheses 

Definition 2.1 A family {Ai}™^ of vector fields in R" with Ai = J2^=i ^ij^j 
called non-degenerate if there exist constants N and C such that for every x G R" and 
every vector v G R" one has the bound 

m 

\\vr<c{i+\\xrfT.('^dx),vr, 

2=1 

with {Ai(x), v) = I]"=i Aj(a;) Vj. 

The conditions on K which we will use below are taken from the following list. 

a. The vector fields Xj with j ~ 0, . . . , m belong to Polf^ and the function / 
belongs to Pol^. 

5o. There exists a finite number M such that the family consisting of {^i}™0' 
{{Xi, Xj]}™ ^Q, { [[Xi, Xj], Xfe] }^ . and so on up to commutators of rank 
M is non-degenerate. 

hi. There exists a finite number AI such that the family consisting of {X^}™ j^, 
{[X;, Xj]}™-^Q, { [[Xi, Xj], Xfe] } . . and so on up to commutators of rank 
M is non-degenerate. 

The difference between 6o and hi is in the inclusion of the vector field Xq (in ho), so 
that hi is stronger than b^. 

Definition 2.2 We call JCq the class of operators of the fonn of (1.1) satisfying a and 
ho above, and tCi the class of those satisfying a and hi. Clearly, hi is more restrictive 
than ho and therefore K-i C /Co- 

Remark 2.3 If K is in /Co then K is hypoelliptic. If K is in /Ci then dt K is 
hypoelliptic. 
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3 Localized Bound 

The main result of this section is Theorem 3.1 which provides bounds for localized test 
functions. 

We let B{x) denote the unit cube around x e R" : 

B{x) ={ye R" I \yj-Xj\<l,j^l,...,n} . 

To formulate our bounds, we introduce the operator A, defined as the positive square 
root of = 1 — X]r=i ^ ^ ^ Later on, we will also need the multiplication 
operator A defined as the positive root of (multiplication by) A^ = 1 + 

Theorem 3.1 Assume K G /Ci. Then, there exist positive constants e*, C*, and N<, 
such that for every x e R" and every u € {B{x)), one has uniformly /or y e R".- 

\\K''u\\<C,{l + \\xff'\\u\\ + \\{K + iy)u\\ . (3.1) 

If K is in /Co (but not in JCi) the same estimate holds, but the constant C* will depend 
generally on y. ' 

Proof. The novelty of the bound is in allowing for polynomial growth of the coeffi- 
cients of the differential operators. Were it not for this, the result would be a special 
case of Hormander's proof of hypoellipticity of second-order partial differential opera- 
tors [Hor85, Thm. 22.2. 1]. Since the coefficients of our differential operators can grow 
polynomially we need to work with weighted spaces. 

We introduce a family of weighted Sobolev spaces 3°"'^ with a,j3 G R as the 
following subset of tempered distributions S,' on R": 

S'"''' = {u e S,; I K^'K^u e l2(R")} . 

We equip this space with the scalar product 

(/,g)„,^ = (A"A''/,A"A^'5)L. , (3.2) 

writing also (•, •)„ instead of (•, •)q.o. We also use the corresponding norms |j • ||a,/3- 
Note that these spaces are actually a particular case of the more general class of Sobolev 
spaces introduced in [BC94]. 

The following lemma lists a few properties of the spaces 3°"'^ that will be useful in 
the sequel. We postpone its proof to Appendix A. 

Lemma 3.2 Let a, /3 G R. We have the following: 

a. Embedding: For a' > a and fj' > f3, the space S" is continuously embedded 
into 3"'^. The embedding is compact if and only if both inequalities are strict. 

b. Scales of spaces: The operators A"^ and are bounded from 3°"'^ into 3°"^^'^ 
and 3°''^~'^ respectively. If X ^ Pol^ then X is bounded from 3"'^ into 

ga~k,l3-N 



The norms are 1? norms. 
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c. Polarization: For every 7, (5 G R, one has the bound 

{f,9)o..0<C \\f\U'.s3' ll.9ll«",/3" , a' + a"^2a, /?' + /?" = 2/? , 

which holds for all f and g belonging to the Schwartz space Sn. The constant C 
may depend on the indices. 

d. Commutator: Let X £ Polf and Y e Pol^'. For every 7 e R, [X, AT'] is 
bounded from S'^'P into ^-a+i-fe-T^/^-w. Similarly [X,[Y,K'<]\ is bounded 
from S'^'P into 5'a+2-fe-fe'-7,/3-w-w'_ 

e. Adjoint: Let X G Pol^ and let f^g^Sn- Then 

{f,Xg)^^p = {X^f,g)^^p + R{f,g), 
where the bilinear form R satisfies the bound 

\R{f,9)\<C\\!\U',3'\\g\U".p" , 

with 

a' + a" ^2a + k-l, fi' + fi" ^ 2(i + N . (3.3) 
The constant C may depend on the indices. 

Notation 3.3 We write Ky instead of K + iy. We also introduce the notation $ < B 
to mean: There exist constants C and N independent of x and y such that for all 

•fr^cii + MfiM + WKyuW). 

We will show below that 

WAA^-^uW < B , (3.4) 

holds for A taking values among all of the vector fields appearing in 61 or 69. Assuming 
(3.4) one completes the proof of Theorem 3.1 as follows: Notice that if the collection 
is non-degenerate, then 

k 

for every x S R" and every u £ C^{13{x)). Therefore, by (3.4) we find 

\\u\\l = \\AA^-'u\\ < Ci(l + ||x|p)^^^||A,A-M' < B^. 

1=1 

Polarizing, we obtain: 

h|le'/2 < MM. < C2Mil + MY^M + \\KyU\\) 

< Ci\\uf{l + \\xff''^- + i\\u\\ + \\Kyu\\f 

< (C2||^.||(l + i|xi|2)^^ + ||yi| + ||A>||)', 

and hence (3.1) follows with = e/2, = N2, and C* = C2 + 1. 
It remains to prove (3.4). 
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Remark 3.4 To the end of this proof, we use the symbols C and N to denote generic 
positive constants which may change from one inequaUty to the next. 

By the bound on [A,K^^^] of Lemma 3.2(d) — and the fact that u e C^{B{x)) 

implies ||u||o,Ar < C{l + for every N > — we will have shown (3.4) 

if we can prove 

\\Au\U-i<B. (3.5) 
Notice that by Lemma 3.2(b), the estimate (3.5) yields 

\\Au\\U^^ < C,(l + \\x\\y+''{\\ur + \\Kyu\\') , (3.6) 

for every 7 > 0, a; G R", and u G C^{B(x)). 

To prove (3.5), we proceed as follows. First, we verify it for A = Xi with i = 
1, . . . , m (as well as for A — Xq in the case /Co). The remaining bounds are shown by 
induction. The induction step consists in proving that if (3.5) holds for some A G Polf' 
then 

||[A^,]"lle/8-i <B for z = 0,...,m. (3.7) 

The first step. By the definition of K and the fact that Xi maps C^{B{x)) into itself, 
we see that 

||X,u||<B, i = l,...,m, (3.8) 

that is, (3.5) holds for e < 1 and A = Xi. 

We next show that it also holds for A = Xq whenever e < 1/2. (This will be the 
only place in the proof where C depends on y, but we need this estimate only for the 
case ICq.) Using (1.1) and Lemma 3.2(c), we can write 

m 

\\X0u\ty2 < \\Xou\\-i{\\Kyu\\ + \\fu\\ + \y\ \\u\\) + Y,{Xou,XjX^u)^^/^ . 

i=l 

Using Lemma 3.2(b) to estimate ||Xo?i||_i, the first term is bounded by B^, so it re- 
mains to bound {Xqu, XJ X^u) -1/2- Using this time Lemma 3.2(e), (with a — — ^ 
and (3 = 0), we write 

{Xou,XjX^u)^y2 = {X,Xou,X,u)_y2 + R{Xou,X,u) , (3.9) 

where R(Xou, Xiu) is bounded by C||Xou||-ij|XiM||, which in turn is bounded by 
B^, using the previous bounds on j|XoMj|-i and The first term of (3.9) can be 

written as 

\{X,Xou,X,u)_,^2\ < C\\X,Xou\\_,\\XM\ ■ 

Since < B by (3.8), we only need to bound ||XjXoM||_i by B. This is achieved 

by writing 

\\X,Xou\\-i < \\XoXM\-i + \\[X,,Xo]u\\^i . 

The second term is bounded by B using Lemma 3.2(b). The first term is also bounded 
by B since \\Xiu\\o,N < C(l + m|| and Xq is bounded from 5*°^^ into 

5"^'° (for some A^) by Lemma 3.2(b). Therefore, we conclude that 

||^oii||-i/2 <B, (3.10) 
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where C will in general depend on y. 

The inductive step. Let A G Polf and assume that (3.5) holds. We show that a similar 
estimate (with different values for e, C, and A^) then also holds for B = [A, X^] with 
i = 0, . . . ,m. We distinguish the case i = from the others. 

The case i > 0. We assume that (3.5) holds and we estimate for some 

e' < 1/2 to be fixed later. We obtain 

\\Buf^,_, = {Bu, AX, u),,^i - {Bu, X,Au)e'-i = Ti + ■ 

Both terms Ti and T2 are estimated separately. For Ti, we get from Lemma 3.2(e): 

Ti = -{ABu,X,u)^'^i + R{Bu,Xiu) , 

where (since e' < 1/2), 

\R{Bu,Xiu)\ < (17(1 + i|a;||)^||Bu||_i||X,u|| < C(l + u|| < . 

(3.11) 

The term {ABu, Xi u)e'-i is written as 

\{ABu,X,u),,.i\ < \\BAu\\2,'^2\\X,u\\ + \\[A, B]u\\^i\\X,u\\ . 

The second term is bounded by like in (3.11). The first term is also bounded by B^ 
by combining Lemma 3.2(b) with the induction assumption in its form (3.6) (taking 
2e' < e). The estimation of T2 is very similar: we write again 

T2 = -{X,Bu, + R{Bu, Au) . (3.12) 

The first term is bounded by C||Xii3u|| _i|j^u||2e'_i. The second factor of this quan- 
tity is bounded by B by the inductive assumption, while the first factor is bounded 
by 

\\X,Bu\\^i < \\BX,u\\^i + \\[B,X,]u\\-i < B , (3.13) 

using Lemma 3.2(b) and the estimate j|Xj||o,jv < B. The remainder R of (3.12) is 
bounded by 

\RiBu,Au)\ < \\Bu\\^i\\Au\\2e'-i , 

which is bounded by B^, using Lemma 3.2(b) for the first factor and the inductive 
assumption for the second. Combining the estimates on Ti and T2 we get 

\\Bu\\,,^i<B for e'<£/2, 

which is the required estimate. 

The case i = 0. To conclude the proof of Theorem 3.1, it remains to bound ||-Bu||£'_i 
by B. In this expression, B = [A, Xq] and e' > is to be fixed later. We first introduce 
the operator 

ni 

k = Y^xlx,, 

i=l 
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which is (up to a term of multipHcation by a function) equal to the real part of Ky, 
when considered as an operator on ]3. We can thus write Xq as 

Xo^K-K + fi=K- + /2 , 

for two functions /i, /2 G Pol^ for some A^. This allows us to express B as 

B = [A, Xo] = AKy + A'jA + [K, A] - 2k A + Afi - /s^ . 

We write j|-Buj|^/_]^ = {Bu, [A, Xo]M)e'_i and we bound separately by each of 
the terms that appear in this expression according to the above decomposition of the 
commutator. 

The two terms containing /i and /2 are bounded by using the inductive assump- 
tion. We therefore concentrate on the four remaining terms. 
The itrm AKy. We write this term as 

{Bu,AKyu)^^^i = -{BAu,Kyu)s>^i + {[A, B]u, Kyu)^'-! + R(Bu, Kyu) , 

where the two last terms are bounded by B^ using Lemma 3.2(b,e). Using assumption 
(3.6) (assuming e' < e/2) and Lemma 3.2(b,c), we also bound the first term by B^. 
The term KjA. We write this term as 

{Bu,KjA),,^i = {KyBu,A),,^i + {A^-^^' [K, A^^' -^]Bu, Au),,^i = Ti + . 

The term Ti is bounded by ||A'j,i3u||_i||y4?i||2e'-i by polarization. The second factor 
of this product is bounded by B, using the induction hypothesis and the assumption 
e' < e/2. The first factor is bounded by 

\\KyBu\\^i < j|BA>j|_i + \\[K,B]u\\_i . (3.14) 

The first term of this sum is obviously bounded by B. The second term is expanded 
using the explicit form of K as given in (LI). The only "dangerous" terms appearing 
in this expansion are those of the form || [X^X^, B]u\\^i. They are bounded by 

\\[Xj'X^,B]u\\^,<\\[Xj',B]XM\-i + \\[X^,B]Xj'u\\^, + \\[X^AX^,B]]u\\^^ . 

The terms in this sum are bounded individually by B, using the estimates on \\Xi u\\, 
together with Lemma 3.2(b,d). We now turn to the term T2. We bound it by 

1^2! < C\\A^-^''[K,A^''-^]Bu\\^i\\Au\\2e'-i . 

The second factor is bounded by B by the induction hypothesis, so we focus on the first 
factor. We again write explicitly as in (1.1) and estimate each term separately. The 
two terms containing Xq and / are easily bounded by B using Lemma 3.2(b,d). We 
also write Xf X^ = Xf + Yi with Yi £ Pol^ and similarly bound by B the terms in 
Yi . The remaining terms are of the type 

= \\A^-^''[Xf,A^''-^]Bu\\-i . 
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They are bounded by 

< 2||A2-2^'[X„A2^'-2]X,B«|1_i + ||A2-2^' [X,,[X^,A^''-^]]Bu\\_^ . 

In order to bound the first term, one writes XiB ^ BXi + {Xi, B] and bounds each 
term separately by B, using the bound ||Xiu||o,-y < B together with Lemma 3.2(b,d). 
The last term is also bounded by B using Lemma 3.2(d). 

The term [A', A]. We write K = X]i=i -^I -^i ^^'^ bound each term separately: 

{BuAX];X,,A\u),,^^ ^ {Bu,X];[X,,A]u),,_^ + {Bu,[Xf ,A]X,u),._^ 
= Tis + Ti^2 ■ 

The first term is written as 

T,,i = {X,Bu, [X^, A\u),,^i + R{u) , 

where R{u) is bounded by C 1 1 Bu 1 1 _ 1 1 1 [X; , A] u 1 1 2e' - 1 . The first factor is bounded by 
B using Lemma 3.2(b) and the second factor is bounded by B, using the estimate for 
the case i ^ (we have to assume e' < e/4 in order to get this bound). The term 
{XiBu, [Xj, is estimated by 

\{X,BuAX^,A]u),,^,\ < \\X,Bu\\^i\\[X^,A]u\\2e'-i ■ 

The first factor is bounded by B as in (3.13) and the second factor is again bounded by 
B, using the estimate for the case i 7^ 0. It thus remains to bound T'i.2, which we write 
as 

T,2 = (Su,X,[Xf , + {Bu, [[Xf ,A],X,]h)^,_^ . 

The first term in this equation is similar to the term {Bu, Xf[X^, A]u)^r^i and is 
bounded by B'^ in the same way. The second term is bounded by 

{Bu, [[Xl,A],X,]u)^,_^ < \\Bu\\^4[[Xl,A],X,]u\\^^,_^ , 

which can also be bounded by B^, using the estimate for the case i ^ 0, provided 

e' < e/8. 

The term KA. In order to bound this term, we need the following preliminary lemma: 

Lemma 3.5 Let v G Sn, a,5 e R, and let Ky be as above. There exist constants C 
and N independent ofy such that the estimate 

m m 

Rc{KyV, v)^ - J2 ll^^^ll" I ^ ^J2\\^'''^L-s,nM^+s,n + C\\v\\l^ , (3.15) 

i=l 1=1 

holds. 

Proof. Obviously (ifyW, w)q = {Kv,v)a. We decompose A' according to (1.1). The 
terms containing Xq and / are bounded by C||wj|^ according to Lemma 3.2(b,e), so 
we focus on the terms containing Xf X^. Using Lemma 3.2(e), we write them as 

{XjX^v,v)^ = \\X,v\\l + R,{v), 

where Ri{v) is bounded by C||Xii)||ct„5_jv|K'||a+5,JV- This concludes the proof of 
Lemma 3.5. □ 
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We now write the term containing KA as 

m 

{Bu,KAu)e'-i =^((X,Bu,XiAu)e'_i + i?i), (3.16) 

i=l 

and we apply Lemma 3.2(e) with / = Bu, g = XiAu, X = Xj. Then we find 

< \\Bu\\^i^N\\X.,Au\\2e'-i < \\Bu\Wj, + \\X,Au\\1^,_, . 

By Lemma 3.2(b), the first term is bounded by B^. Using Lemma 3.2(c) to polarize the 
scalar product in (3.16) we thus get 

m m 

\{Bu,KAu),,^,\ < b2 + \\X:Bu\W + cY,U^M\l'-l ■ 

i=l i=l 

The term involving ||Xi£?u||'L]^ is boundedby as in (3.13). The last term is bounded 
by Lemma 3.5, yielding 

m 

\{Bu,KAu),,^i\ < b2 + C\{KyAu,Au)2e'-i\ +CY,\\X^Au\\\^ 



N 

i=l 



C\\Auf 



4e'-l,7V ■ 

The last term in this expression is bounded by B^ by the induction hypothesis if we 
choose e' < e/4. The term containing XiAu can be bounded by B^ as in (3.13), so 
the only term that remains to be bounded is \{KyAu, Au)2e'-i\- By polarizing the 
estimate obtained by Lemma 3.2(c), one gets 

\{KyAu,Au)2e'-i\ < C\\Au\\l,_^ + C\\KyAu\W . 

The first term is bounded by B^ using the induction assumption. The second term is 
bounded by B^ exactly like (3.14) above. Summing all these bounds this proves (3.7) 
and hence the inductive step is completed. 

Since K was assumed to satisfy ICi (or /Co), we see that after M inductive steps 
the proof of Theorem 3. 1 is complete. □ 

4 Global Estimate 

The results of the previous section were restricted to functions u with well-localized 
compact support. In this section, we are interested in getting bounds for every u E Sn. 
The main estimate of this section is given by 

Theorem 4.1 Assutne K is in JCi or in JCq and let Ky = K + iy be as above. For 
every e > 0, there exist constants S > and C > such that for the norms defined by 
(3.2) one has 

\\u\\s,5<C{\\u\\a,, + \\Kyu\\) (4.1) 
holds for every u € S^. The constants C and 5 are independent of y if K € K-i. 
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Since 5''^ '' is compactly embedded into L^, this result implies: 
Corollary 4.2 Let K be as above. If there exist constants e, C > such that 

\\uh^,<C{\\u\\ + \\Ku\\), (4.2) 
then K has compact resolvent when considered as an operator acting on \3. 
Proof of the Corollary. Combining (4.1) with (4.2), we get 

\\u\\s,s<C{\\u\\ + \\Ku\\) . 

This implies that for A outside of the spectrum of K, the operator (/v — A)^^ is bounded 
from \3 into S^'^ . By Lemma 3.2(a), it is therefore compact. □ 

Proof of Theorem 4.1 . Let and N^, be the values of the constants obtained in esti- 
mate (3.1) of Theorem 3.1. Observe that Theorem 3.1 also holds for any bigger value 
of iV*, and we will assume N^, is sufficiently large. 

We choose e > 0. As a first step, we will show that there exist constants 5 and C 
such that, for any x G R" and u £ C^(B{x)), the following estimate holds: 

II'^IIm- < cil + \\u\u, + c(i + My^^\\u\\ . (4.3) 

Denote by J the smallest integer for which 

J>1 + ^, 

e 

and define 

5 = min{2A^„|,^} . (4.4) 

First, we note that when ^ is a positive self-adjoint operator on some Hilbert space Ti, 
one has the estimate 

\\Auf < C\\A''u\\ \\uf-^ , (4.5) 

whenever both expressions make sense. In the case J = 2^ for j an integer, this can 
be seen by a repeated application of the Cauchy-Schwarz inequality. It was shown in 
[KS59] to hold in the general case as well. 
We next use Jensen's inequality to write 

(1 + !i.in-*+^/^iiA^ii < c(^) Vii + c{i + ii.in'-*+^/^'(i+^)||,„ . 

Dividing this expression 

by + and using the definition of J, we get 

(1 + ii.in^/^iiA^ii < c(i + bin-^* (M^) '\\u\\ 
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Using (4.5), the fact that gf- < by (4.4), and u £ C^(B{x)), we get (4.3). 

In order to prove Theorem 4.1, we use the following partition of unity. Let xo ■ 
R [0, 1] be aC°° function with support in |a;| < 1 and satisfying J^iez Xo(^~^) — 1 
for all a: e R. The family of functions 

V = {x.:-R"^ [0, 1] I X € Z"} , 

defined by 

n 

is therefore a partition of unity for R". By construction, when x,x' E Z then Xx and 
Xx' have disjoint support if there exists at least one index j with \xj — x'j\ > 2. We 
can therefore split V into subsets Vk\k^i 3^ such that any two different functions 
belonging to the same Vk have disjoint supports. 

Consider next an arbitrary function u E Sn- We define Ux = XxU, and then the 
construction of the Vk implies 

J2 ll^-llo.- ^ 3"||^i||o.e . (4.6) 

Using (4.3), then Theorem 3.1 and (4.6), we find 

\\uUs < J2 <cJ2{(^ + MT""' 11^-11- + (1 + M'f'WuxW) 
<cJ2 {\M + (1 + ll^ll')^"^* W^y^^W + (1 + ll^lPr^'ll".!!) 

xGZ" 

< CT(\\u\\ + ll^llo..) + C ^ (1 + )-'^*||A>,|| . 

xeZ" 

For k <E {1, . . . , 3"} we now define 

(i + ikir)''^*xM- 

With this notation, we have 

3" 

\\uUs<C\\u\\o,e + cY,\\Kvhu\\ . 

fc=l 

The claim (4.1) thus follows if we can show that 

\\Kyfku\\ < C\\u\\+C\\Kyu\\ . (4.7) 

Since the fk are bounded functions, it suffices to estimate || [K, fk]u\\. By construction, 
every derivative of fk decays like (l + ||a;|P) 

Note that for sufficiently large N^, the functions [Xj, fk] and [Xk, [Xj, fk]] are 
bounded. Since Theorem 3.1 allows us to choose as large as we wish, (4.7) follows 
from the estimate < ||i^.yu||. □ 
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4.1 Cusp 

Our statement about the cusp-like shape of the spectram of K is now a consequence of 
Theorem 4. 1 . 

Theorem 4.3 Let K £ Pol^ be of the type (1.1). Assume that the closure of K in \3 is 
m-accretive and that K g /Ci. Assume furthermore that there exist constants e, C > 
such that 

\\uh^,<C{\\u\\ + \\Kyu\\), (4.8) 
for all y G R. Then, the spectrum of K (as an operator on \3) is contained in the cusp 

{A e C I ReA > , |ImA| < C(l + ReA)''} , 

for some positive constants C and v. 

Remark 4.4 In principle, our proofs give a constructive upper bound on v. However, 
no attempt has been made to optimize this bound. 

Proof. The proof follows very closely that of Theorem 4. 1 in [HN02] , however we give 
the details for completeness. One ingredient we need is the following lemma: 

Lemma 4.5 Let A : I? ^ I? be a maximal accretive operator that has as a core. 
Assume there exist constants C,a > Ofor which 

\\Au\\ < C\\u\\a,a , Vu e Sn, . 

Then, for every N CzN, there exists a constant Cn such that 
Proof. By Lemma 3.2(b), one can bound jluHa.a by 

The generalized Heinz inequality presented in [Kat61] then yields 

This concludes the proof of Lemma 4.5. □ 

We now turn to the proof of Theorem 4.3. Since K e Pol^, one has for a = 
max{2, N} the bound 

\\{K + l)u\\ < C\\u\U^ , yueSn. 

By Lemma 4.5, one can find for every (5 > an integer A/ > and a constant C such 
that: 

{u, {{K + iTiK + l))'/^'u) < C\\u\\ls , (4.9) 
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Furthermore, Theorem 4.1 together with (4.8) yields constants C and 5 such that for 
every u £ Sn and every y G R: 

\\u\\ls<C{\\uf + \\(K + ty)u\\^) . (4.10) 

Since K is ;7z-accretive by assumption, we can apply [HN02, Prop. B.l] to get the 
estimate 

i|z + Ip/^^hf < {{(K + mK + l)fu,u) + \\(K - z)uf 
< C\\u\\ls + \\iK-z)uf, 

where the second line is a consequence of (4.9). Using (4. 10) and the triangle inequality 

for z = Rez + i Imz, we get 



1 



-\z + i|w||- < C{{l + Rezf\\u\f + \\(K~z)u\f) . 

Together with the compactness of the resolvent of K, this immediately implies that 
every A in the spectrum of K satisfies the inequality 

i|A + l|"/^^||u!P < C(l + ReA)2||u||2 . 

This concludes the proof of Theorem 4.3. □ 

5 Examples 

We present two examples in this section: A first, very simple one, and a second which 
was the main motivation for this paper 

5.1 Langevin equation for a simple anharmonic oscillator 

Our first example consists of one anharmonic oscillator which is in contact with a 
stochastic heat bath at temperature T. The Hamiltonian of the oscillator is given by 

For this model the associated spectral problem can be solved explicitly when e = 0, 
because it is an harmonic oscillator The spectrum lies in a cone as shown in Fig. 2. 
We also show that in first order perturbation theory in e, the spectrum seems to form a 
non-trivial cusp, but this result remains conjectural, because of non-uniformity of our 
bounds. 

The Langevin equation for this system is 



dp = -v'^q dt - eq^ dt -'-fpdt+ y/2-iTdw{t) , dq=pdt, (5.1) 

where 7 > measures the strength of the interaction between the oscillator and the 
bath. Denote by (O, P) the probability space on which the Wiener process wit) is de- 
fined. We write ift.Miix) with lo £ Q, for the solution at time t for (5.1) with initial 
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condition x = (p, q) and realization u) of the white noise. The corresponding semi- 
groups acting on observables and on measures on are given by 

{Ttf)(x)^ f {foipt,Ux))dF(Lo), (5.2a) 

{T;y){A)= f {fioip-^(A))dP(u;), (5.2b) 
Jn 

where A C is a Borel set. It is well-known that 

dfirip, q) = exp {-H(p, q)/T) dpdq 

is the only stationary solution for (5.2b). 

The Ito formula yields for ft ^ Ttf the Fokker-Planck equation given by 

dtft = iTdlft +pdjt- + eq' + 7P) dp ft . (5.3) 

We study (5.3) in the space Tip ~ \3(^ ^df^ir)- and make the change of variables 
ft = exp {H/(2T))Ft in order to work in the unweighted space Tio — L^(R^, dpdq). 
Equation (5.3) then becomes dtFt = —C^Ft, where the differential operator is 
given by 

4 = -iTdl + -^p' -l-pd^ + u\dp + eq' dp . 
By rescaling time, p and q, one can bring to the form 

Ce ^ ^{-dp + p^ - 1) + a{qdp- p dq) + ceq' dp , 

where a = 2\pxTv and c > 0. 

The operator K = is thus of the type (1.1) with Xo = a{q dp — p dq) + ceq'^ dp 
and Xi = dp. We now verify the conditions of Section 2.2. It is obvious that 
these vector fields are of polynomial growth, thus condition a is satisfied. Since 
[Xi,Xo] = —adq, the operator satisfies condition bi as well, and so the con- 
clusion of Theorem 4.1 holds. Proceeding like in [EHOO, Prop. 3.7], one shows an 
estimate of the type (4.8) (see also the proof of Theorem 5.5 below, where details are 
given). Therefore, Theorem 4.3 applies, showing that the spectrum of is located in 
a cusp-shaped region. In fact, we show in the next subsection that the cusp is a cone 
when e = 0, and then we study its perturbation to first order in e. 



5.1.1 First-order approximation of the spectrum of £e 

We will exphcitly compute the spectrum and the corresponding eigenfunctions for Cq 
and then (formally) apply first-order perturbation theory to get an approximation to the 
spectrum of C^. We introduce the "creation and annihilation" operators 

P + dp , p-dp q + dq q-dq 
a = , a = ^ , — ^ — 
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in terms of which can be written as 

£e = a* a + Q!(6*a — a*h) + ceq'^ dp . 

With this notation, it is fairly easy to construct the spectrum of £o- Note first that is 
an eigenvalue for £ with eigenfunction exp(— — This is actually the vac- 

uum state for the two-dimensional harmonic oscillator in quantum mechanics (which 
is given by a*a + 6*6), so we call this eigenfunction 

A straightforward calculation shows that the creation operators Cj_ defined by 



4 = a* +(3±h* , 
satisfy the following commutation relation with Cq 



1 



V4a2 - 1 
2a 



1 \/Aa^ - 1 
[£o,cy = A±c^, A± = -±i 



a 



Therefore, Aq = nA+ + mA_ with n and m positive integers are eigenvalues for Cq 
with eigenvectors given by 

{cxr{c*_r\n) . 

We conclude that for a > 1/2 the spectrum of Cq consists of a triangular grid located 
inside a cone (see Figure 2). 



ImA 



ImA 



ReA 



■ReA 



Figure 2 Spectrum of Co. 



Figure 3 Approximate spectrum of Ce 



Remark 5.1 Although the spectrum of £o is located inside a sector, £o is «of sectorial 
since the closure of its numerical range is the half -plane ReA > 0. 

In order to do first-order perturbation theory for the spectrum of we also need 
the eigenvectors for Cq, which can be obtained by applying successively d*^ and d*_ to 
|ri), where 



d*^= a* — P^b* 



With this notation, (d!j_)"(fil)™|0) is an eigenvector of Cq with eigenvalue Aq' 
first-order perturbation theory, the eigenvalues of are approximated by 



A'J 



A" 



(f7|d"d"g^9p(cp"(cl)"|fl) 
(17|<i!"d'|(c;)"(cl)"|17) 



By 



(5.4) 
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The resulting spectrum^ is shown in Figure 3 (the sector containing the spectrum of Cq 
is shown in light gray for comparison). One clearly sees that the boundary of the sector 
bends to a cusp. A (lengthy) explicit computation also shows that 

^„.o = -12n(n - 1) , ^ + 9n . 

In principle this confirms the cusp-like shape of the boundary, were it not for the non- 
uniformity of the perturbation theory (in n). 

5.2 A model of heat conduction 

In this subsection, we apply our results to the physically more interesting case of a 
chain of nearest-neighbor interacting anharmonic oscillators coupled to two heat baths 
at different temperatures. We model the chain by the deterministic Hamiltonian system 
given by 

N 2 N 

(We will give conditions on the potentials Vi and V2 later on.) In order too keep nota- 
tions short, we assume pi, qi e R, but one could also take them in R'^ instead. The two 
heat baths are modeled by classical free field theories and with initial condi- 
tions drawn randomly according to Gibbs measures at respective inverse temperatures 
(3l and Pr. (We refer to [EPR99a] for a more detailed description of the model.) It is 
shown in [EPR99a] that this model is equivalent to the following system of stochastic 
differential equations: 

dqi = pidt , i = 0, . . . , N , 

dpo = -Vliqo) dt + V^(qi) dt + rLdt, 

dp, = -Vliq,) dt - V^(q,) dt + V^(qj+i) dt , j = 1, . . . , TV - 1 , 

dpN = ~V{{qN) dt - V^iqN) dt + rB,dt , 



drL = -Jltl dt + XijLqo dt - Xl\/'2jlTl dwL(t) , 



dm = -7i?Tfl dt + X%"fRqN dt - \R^j2^iiTRdwR{t) , 

where Ti = /3j^, 7i are positive constants describing the coupling of the chain to the 
heat baths, and Wi are two independent Wiener processes. The variables ?! and 
describe the internal state of the heat baths. If Tl — Tr = T, the equilibrium measure 
for this system is d^xip, f) = exp {—G{p, q, r) /T) dp dq dr, where the "energy" G 
is given by the expression 

G{p, q, r) = H{p, q) + - qo?X + - iNrR ■ 



-Actually the set {A^ + ceSn,m \ n,m > 0}. 
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If Tl ^ Tn, there is no way of guessing the invariant measure for the system. We can 
nevertheless make the construction of Section 5.1 with the reference measure dfif for 
some temperature 

r>max{Ti,T^}, 

which is a stability condition, as one can see in (5.6) below. The resulting operator 
= £ is given by 

C = XIX^ + X*j,Xj, + fl + fl + Xo, (5.5) 

where 



fh.R = V 1l.r{Tl,rIT ~ l){rL,R - Xl.rQo.n) , (5.6) 



Po 



with 



Pn 



'^L,R^-^^(TL,R^f) 
■^L,R 



We are now in a position to express the conditions of Section 2.2 in terms of sufficient 
conditions on the potentials of the model. The first two assumptions guarantee that C 
is in /Ci. 

Assumption 1 There exist real numbers n,m > such that D"Vi <E Polg"^" and 
D°'V2 e Polg'""" /or a < 2. 

Assumption 2 There exists a constant c > such that V2'{x) > cfor all a; e R. 

Remark 5.2 The second assumption states that there is a non-vanishing coupling be- 
tween neighboring particles in every possible state of the chain. 

The verification that these assumptions imply a is easy, and the verification that bi 
holds can be found in [EPR99a, EHOO]. 

Proposition 5.3 Let C be defined as above and let Vi and V2 fulfill Assumptions 1 and 
2 above. Then C satisfies the assumptions of Theorem 4.1 and satisfies Eq.(4.1) with C 
and 5 independent ofy. 

In order to show that the spectrum of C is located in a cusp-shaped region {i.e. that 
the hypotheses of Theorem 4.3 hold), two more assumptions have to be made on the 
asymptotic behaviour of Vi and V2 : 



Assumption 3 The exponents n and m appearing in Assumption 1 satisfy 1 < n < m. 
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Remark 5.4 The physical interpretation of the condition n < m (actually 1 < n < m 
would probably work as well, see [RBT02b], but we could not apply directly the results 
of [EHOO]) goes as follows. If n > m, the relative strength of the coupling between 
neighboring particles decreases as the energy of the chain tends to infinity. Therefore, 
an initial condition where all the energy of the chain is concentrated into one single 
oscillator is "metastable" in the sense that the energy gets transmitted only very slowly 
to the neighboring particles and eventually to the heat baths. As a consequence, it is 
likely that the convergence to a stationary state is no longer exponential in this case, 
and so the operator C has probably not a compact resolvent anymore. 

Our last assumption states that the potentials and the resulting forces really grow asym- 
ptotically like \x\^^ and respectively (and not just "slower than"). 

Assumption 4 The potentials Vi and V2 satisfy the conditions 

Vl{x) > Ci(l + WxWy - C2 , xVlix) > C3(l + WxfY - C4 , 

V2(X) > C5(l + \\x\fr - C6 , ^V^(X) > C7(l + - C8 , 

for all X €z R and for some positive constants Ci. 

Theorem 5.5 Let C be defined as above and let Vi and V2 fulfill assumptions 1^ 
above. Then, C has compact resolvent and there exist positive constants C and N such 
that the spectrum of C is contained in the cusp 

|a e C ReA>0 and ImA < C(l + |ReA|)"| . 

Proof. We will apply Theorem 4.3, and need to check its assumptions. It has been 
shown in [EHOO, Prop. B.3] that C is m-accretive. The fact that C £ K.i was checked 
above, and (4.8) was shown for ?/ = in [EHOO, Prop. 3.7]. However, closer inspection 
of that proof reveals that whenever was used, it only appeared inside a commutator. 
Therefore, we can replace it by Xq + iy without changing the bounds. Thus, we 
have checked all the assumptions of Theorem 4.3 and the proof of Theorem 5.5 is 
complete. □ 



A Proof of Lemma 3.2 

The points a and b of Lemma 3.2 are standard results in the theory of pseudodifferential 
operators (see e.g. [Hor85, Vol. Ill] or, more specifically, [BC94, HT94a, HT94b]). The 
point c is an immediate consequence of the Cauchy-Schwarz inequality combined with 
a. In order to prove the points d and e, we first show the following intermediate result: 

Lemma A.l Let f : R" R and a G R. Let k be the smallest even integer such that 
I a I < k. Then, if f satisfies 

sup \d^f{y)\ <K, y \S\<k, 

the corresponding operator of multiplication is bounded from S"'^ into S"'^ and its 
operator norm is bounded by Ck. The constant C depends only on a and (3. 
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Proof. By the definition of 3°"'^ , it suffices to show that the operator A"fA~°' is 
bounded by Ck from into L^. Since / is obviously bounded by k as a muhipU- 
cation operator from into \J, it actually suffices to bound A"[/, A^"]. Assume first 
that a £ (0, 2). In that case, we write 

poo I. a 1 

A"[/,A-"] = a/ z-^/'—-j^[f,A']—j^dz. 

Jq z + z + A^ 

The commutator appearing in this expression can be written as 

n 

[/,A2]=^(29J9, + 92/) . (A.1) 

i=l 

It is clear from basic Fourier analysis that \\di{z + A^)^^/^ j| < 1 and therefore 

||[/,A2](z + A2)-i/2|| <Ck. 

Furthermore, the spectral theorem tells us that for any function F, || F(A^)|| is bounded 
by sup_,^>]^ F{X). Therefore there exists a constant C independent of z > such that 

||A"(z + A2)-1|| < 



Combining these estimates shows the claim when a e (0, 2). The case a ~ 2 follows 
from the boundedness of [/, A^] A"^. Values of a greater than 2 can be obtained by 
iterating the relation 

A"+VA""^^ = A"/A-" + A"[/, A2]A-"-2 

Using (A.l), the fact that di commutes with A, and the fact that 9; A"^ is bounded, we 
can reduce this to the previous case, but with two more derivatives to control. The case 
a < follows by considering adjoints. This concludes the proof of Lemma A.L □ 

Remark A.2 Since the direct and the inverse Fourier transforms both map S'"''' con- 
tinuously into S*^'", the above lemma also holds for bounded functions of dy and not 
only for bounded functions of y. 



We are now ready to turn to the 

If 



Proof of point d. Let X S Polf . We first consider 7 € (—2, 0). Since, in Fourier 
space, A^ is a multiplication operator by a real positive function, we can write 

{X, A^] = / z''/^ ^— [X, A\ 



A2 ' Z + A2 ■ 

In order to bound this expression, we define B = {X, A^], commute B with the resol- 
vent, and obtain 

[X,A'']^cJ z^/2A_^A^-25 + C, / -^——[B,A']-' 



(Z + A2)2 ^ (Z + A2)2^ ' ■'Z + A2 
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The first term equals C:I|,A''' because z'^/^a;^ '^(z + x^) ^ does not depend 
on a; > 0. This, in turn, is bounded from S'"''^ into 5'a+i-fc-7,/3-N ^^jj^g ^ g Polf+i 
and Lemma 3.2(b). To bound the second term, we rewrite 

The factor A^(z + A^)^^ is bounded from 3°^'^ into itself, uniformly in z. Using 
Lemma 3.2(b) as before, we see that the factor K''~^{B^ A^]A^-^ is bounded from 3°"'^ 
into sa+i-k-'i,i3-N ^ 5"" ,/3 Finally, using Lemma A.l and counting powers, we 
see that the first factor has norm bounded by 0(z-3/2) f^j. j^^gg ^ ^j^^j 0(^:7/2) fo^ z 
near as a map from 3°" to itself. This proves the first statement of Lemma 3.2(d). 
The second one is proven similarly and is left to the reader. 
Proof of point e. Recall that we want to bound 

I=\{f,Xg)^.p-{X^f,g)^,0\, 

where X e Pol^ and X'^ denotes the formal adjoint (in I?) of X. We write this as 

I={[A-^^'A-'-,X^]A^-A'Pf,g)^,f,. 

We rewrite the operator as 

[A-2/9A-2",X^]A2"A2'3 = A-2'3[A-2«,X^]A2"A2'3 + [A'^^', X^] A^^ . 

The second term is in Pol^_]^ by inspection, and the required bound follows at once 
from Lemma 3.2(b,c). The first term is bounded similarly by using Lemma 3.2(d,b,c). 
This concludes the proof of Lemma 3.2. □ 
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